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?n ■ ABSTRACT 



The recent detection of gamma radiation from Mkn 501 at energies as high 
as ~ 25 TeV suggests stringent upper bounds on the diffuse, far infrared, 
extragalactic radiation density. The production of electron-positron pairs 
O ! through photon-photon collisions would prevent gamma photons of substantially 

Q>^ ! higher energies from reaching us across distances of order 100 Mpc. However, 

^ ! coherently arriving TeV or sub- TeV gammas — Bose-Einstein condensations of 

^! photons at these energies — could mimic the Cherenkov shower signatures of 

2 ! extremely energetic gammas. To better understand such events, we describe 

to '. their observational traits and discuss how they might be generated. 

> 

^ ■ Subject headings: BL Lacertae objects: individual (Mkn 501) — diffuse radiation 

■ — gamma rays: theory — infrared: general — masers 



1. Introduction 

High energy gamma rays are readily absorbed in the intergalactic medium through 
pair production in a sufficiently dense, diffuse, microwave or infrared radiation field (Gould 
& Schreder, 1966; Stecker, De Jager, & Salamon 1992). For this reason, a great deal of 
attention has be paid to gamma rays at energies apparently reaching ^ 10 TeV, recently 
detected from the galaxy Mkn 501 (Hayashida et al., 1998, Plan et al., 1998, Aharonian 
et al., 1999, Krennrich, et al., 1999). Mkn 501 is a BL Lac object at a distance of 
~ 200 Mpc, for a Hubble constant, Hq = 50 kms~^Mpc~^. Unattenuated transmission 
of ^ 10 TeV photons across distances of this order would place severe constraints on the 
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diffuse extragalactic infrared background radiation (Coppi & Aharonian, 1997, Stanev 
& Franceschini, 1998) placing upper limits to the radiation density that are close to 
values derived from COBE detections and IRAS source counts alone (Hauser, et al., 1998; 
Hacking & Soifer, 1991; Gregorich, et al., 1995). Given these close coincidences it is useful 
to re-examine the severity that these observations place on the density of the diffuse 
extragalactic infrared radiation (DEIR). 

2. Bose-Einstein Condensations of Photons 

Coherent radiation, i.e. highly excited quantum oscillators, are produced in a variety of 
processes, but are also regular components of blackbody radiation in the Rayleigh- Jeans tail 
of the energy distribution. These excited oscillators correspond to densely occupied radiation 
phase cells — a Bose-Einstein condensation of photons all having quantum-mechanically 
indistinguishable properties, i.e. identical momenta, positions, polarizations, and directions 
of propagation, within the Heisenberg uncertainty constraints. 

Given that cosmic ray particles can have energies going up to 3 x 10^'^ cV, and given 
that one expects a cutoff for gammas from Mkn 501 at energies many orders of magnitude 
lower, around 10 or 20TeV, it does not seem far-fetched to think that the actually observed 
gammas reaching Earth might lie far out in the low-frequency tail of some significantly 
more energetic radiation field characterized by an equivalent temperature much higher than 
a few TeV. 

If this were the case, we would expect that the radiation arriving at Earth could be 
highly coherent, meaning that phase cells would be filled to rather high occupation numbers, 
N. As they interact with the DEIR, densely filled phase cells can decline in population and 
lose energy only by going stepwise from an initial occupation number A^, to (A^ ~ 1)? and 
from there to (A^ — 2), etc. Because the mean free path for interactions of photons with the 
DEIR is energy dependent, a fraction of a coherent assembly of photons could penetrate 
appreciably greater distances through the diffuse extragalactic radiation field than, say, a 
single photon of the same total energy. 

A number Na of such arriving photons, each with energy hu would impinge on the 
Earth's atmosphere at precisely the same instant, and would interact with the atmosphere 
producing an air shower that emits Ghcrcnkov light that could mimic that due to a single 
photon with energy Nahv impinging on the atmosphere. These two kinds of impacts could 
be distinguished by shower images they produce and probably also by the fluctuations in 
the energy distribution observed near the cut-off energy Eqq for a series of Cherenkov 



-3- 



events. 

Because of their high momenta, the arriving bunched photons would spread over 
only the smallest distance Ay in their traversal through extragalactic space, given by the 
uncertainty relation ApyAy ~ h, where Apy is the uncertainty in transverse momentum. 
Apy is the product of the photon momentum hu/c and the angular size that the source 
subtends at Earth. The smallest dimension wc could expect would be of the order of 
an AGN black hole Schwarzschild radius ~ 3 x 10^^ M / {10^ Mq) cm. This would make 
Ay ~ (10^Mq/M)10~^ cm — negligible in Cherenkov detection. 

3. Interpretation of Cherenkov Radiation Data 

TeV 7-rays are detected through the Cherenkov radiation generated in the Earth's 
atmosphere by electrons in an "air shower" initiated by the 7-ray. Such air showers are 
electromagnetic cascades involving pair production and bremsstrahlung interactions. As 
long as the energy of the photon entering the atmosphere is sufficiently high, the Cherenkov 
yield of the air shower is sensitive primarily to the total energy deposited, not to the number 
of instantaneously arriving photons. Accordingly, one might expect such telescopes to 
mistakenly record five simultaneously arriving 5 TeV photons as a single shower of 25 TeV. 
On the other hand, if the number of simultaneously arriving photons, N, were much higher, 
then the showers would look very different, and if N were really large there would be no 
Cherenkov radiation at all. 

To quantify the discussion above, we shall compare the mean and standard deviation 
of the number of electrons in the shower, Ne{t), as a function of depth into the atmosphere 
measured in radiation lengths, t, for the two cases. Note that the atmosphere is 
approximately 1030 g cm~^ thick and the radiation length of air including forward 
scattering is 36.66 g cm^^. Although the cross section for interaction of an assembly of 
coherent photons is A^ times higher than that of an individual photon, a shower initiated 
by an assembly of A^ coherent photons having total energy Ne would be identical to a 
superposition of A^ showers due to individual photons of energy e. Above ~ 3 GeV the pair 
production mean free path for photons in air is constant at ipair — 9/7 radiation lengths. 
For an assembly of N coherent photons, the pair production mean free path is therefore 
identical to an exponential distribution with mean tpgXv — (9/7), i.e. it is the same as the 
distribution of first interaction points of a single photon. This also implies that at depth t 
the average number of photons remaining in the assembly is N exp{—t / tpgXv) ■ 

Crewther and Protheroe (1990) provide a parametrization of the distribution of the 
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number of electrons in photon initiated showers, p[Ne{t — ti)], as a function of depth into 
the atmosphere beyond the first interaction points of the primary photons, (t — ti). We 
use their results together with our Monte Carlo simulation of the first interaction points 
of each of the N photons in a coherent assembly to simulate the development of the air 
shower due to the coherent assembly, thus taking account of all fluctuations in shower 
development. In Fig. we show as a function of atmospheric depth t Ng and ± la 
based on 1000 simulations for the case of single photons of energy 25 GeV, assemblies of 5 
coherent photons each having energy 5 TeV, and assemblies of 25 coherent photons each 
having energy 1 TeV (each assembly has energy 25 GeV). As can be seen, air showers 
due to coherent assemblies develop higher in the atmosphere, and have much smaller 
fluctuations in shower development. Such differences between showers due to single photons 
and assemblies of coherent photons would produce different Cherenkov light signatures and 
should be detectable with state-of-the-art Cherenkov telescopes such as HEGRA (see e.g. 
Konopelko et al. 1999). 

4. Extragalactic Optical Depth 

Propagation of assemblies of A^o coherent photons each of energy e through the 
microwave and DEIR fields is analogous to their propagation through the atmosphere. 
However, assemblies of coherent photons having total energy ii^tot = N^e may travel farther 
than single photons of energy N^e without interaction because, unlike in the atmosphere, 
the mean free path for pair-production in the extragalactic radiation fields depends strongly 
on photon energy. 

Just as in the air-shower cascade, only a single photon at a time can be lost from 
a phase cell, with a corresponding decline in occupation number from A^ to A^ — 1. On 
each encounter with an infrared photon, the coherent assembly of A^ photons has an 
A^-fold increase in probability for some photon to be removed, so the mean free path is 
Xpi,ir{e)/N where 

2^pair(^) is the mean free path for photon-photon pair production by 
single photons of energy e through the extragalactic radiation fields. This implies that at 
distance x from the source the average number of photons remaining in the assembly is 
Nr{x) = Noexp[—x/xpa.i^{e)], precisely the expression that would hold for A'o independent 
photons. If d is the distance from the source to Earth, then the energy observable by 
Cherenkov telescopes is £^obs = NFi,{d)e, and the number of photons in the assembly of 
coherent photons on emission was Aq = A'Rexp[(i/a:pair(-E'obs/A^R)]. For the purpose of 
illustration, we use for Xpair(-E') the logarithmic mean of the upper and lower curves of 
Fig. 1(a) of Bednarek and Protheroe (1999) which is based on the infrared background 
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models of Malkan and Stecker (1998). We show in Fig. ^ the result for propagation of 
coherent photons through the microwave and DEIR fields across d = 200 Mpc appropriate 
to Mkn 501. We note, for example, that a coherent assembly of forty 10 TeV photons 
emitted would typically arrive at Earth as a coherent assembly of ten 10 TeV photons with 
an observable energy of 100 TeV, while a single photon of 100 TeV would have a probability 
of much less than 10~^ of reaching Earth. 

5. Fluctuations in the Arriving Phase Cell Energy Content 

A stream of photons characterized by a brightness temperature Tf, of necessity will also 
have a distribution of phase cell occupation numbers, A^, which, for high average values 
{N) fluctuates as (AA^)i.ins ~ (A^). For emission of a stream of identical assemblies of 
coherent photons, each containing Nq photons on emission, fluctuations in the number of 
photons, Nji, remaining in each assembly after propagation to Earth through the DEIR, 
are Poissonian about the mean value (Nr), i.e. (AiV)rms ~ \f{^^^y ^^i {Nji) <^ Nq, and 
less than Poissonian for N^i ~ Nq. Both these effects broaden the energy distributions of 
observed Cherenkov events. 

6. What Mechanisms Could Produce Coherent TeV Gammas? 

In the laboratory (such as DESY), coherent X-radiation can be produced by stimulated 
emission of relativistic electrons through a periodically varying magnetic field (Madey, 
1971). Therefore this shows that such processes are available in principle. 

A more promising astrophysical process might arise from the interaction of a coUimated 
beam of relativistic electrons moving roughly upstream against an OH or H2O megamaser. 
This process is attractive, because a substantial number of AGNs are known to have nuclear 
megamasers. Inverse Compton scattering would produce photons with an energy increase 
7^ in the co-moving frame of the jet of relativistic, randomly directed electrons. Here 7 is 
the Lorentz factor of electrons in the jet's co-moving frame. To produce ITeV photons 
from H2O megamaser radiation at 22 GHz, we would require (57)^ ~ 1.1 x 10^^{E/TeV), 
where 6 = [T{1 — Pcos6)]~^ is the Doppler factor, j3 = v/c refers to the relativistic bulk 
velocity v, and F = [1 — ti^/c^]"^/^ is the Lorentz factor of the jet. The factor 5^ translates 
the photon's initial energy to the co-moving frame and back to the frame of an Earth-based 
observer. For Mkn 501, the line of sight angle 6 appears to be directed very nearly in our 
direction, so we may choose 5 = 2r ~ 25 (e.g. Tavecchio et al. 1998). 
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As shown below, the number of phase cells into which the maser photons can be 
inverse- Compton scattered is limited and quickly fill up for relativistic jets with high column 
densities. At the photon densities discussed, nonlinear effects can be neglected. 

To provide a representative example, we might cite conditions in the galaxy NGC 1052, 
which contains a water megamaser with components that appear to lie along the direction 
of a radio jet (Claussen et al. 1998). Though this may just be a projection effect, we will 
assume as these authors have that it may signify interaction of the jet with dense clumps of 
molecular clouds - possibly producing maser activity in shocks. 

The observed radiation intensity of the maser per unit bandwidth at 22 GHz is 
/(i^) = (cp(i/)/47r) = 50 mJy for a beam size that is unresolved at 0.3 x Imas. The 
beam, however, is clearly much larger than the roughly forty individual sources that are 
detected by virtue of their different velocities along the line of sight, whose centroids are 
separated by as little as ~ 0.1 mas. The brightness temperature of these individual sources 
is Tij{u) = [/(z/)c^/2/cz/^] > 4.5 x 10^ K if the nominal beam size is assumed. The density of 
phase space cells at this frequency is n(z/) = (Sttz/^/c^) ~ 4.5 x lO"^*^ cm~^Hz~^ so that the 
phase cell occupation number becomes A^occ = — (kTb/hu) > 4.3 x 10^. 

All these figures are lower limits, since neither the angular resolution nor the spectral 
resolution suffice to resolve the individual maser sources. For this reason, it may be 
better to assume the properties of the better-resolved Galactic H2O masers, which have a 
brightness temperature of order Th ~ 10^^ K, and a corresponding occupation number of 
order A^occ ~ 10^^ (Moran, 1997). To be somewhat more conservative, we will adopt a value 
of A^occ ~ 3 X 10^3 below. 

Under a Lorentz transformation 7(i/) and p{u) scale as u^, as does hvn{i'), so that the 
phase cell occupation number transforms as a constant. We can therefore deal with the 
occupation number as though it were in the rest frame of the jet of relativistic electrons. 
These electrons with energy ^rrieC^ will have some velocity dispersion, leading to an energy 
bandwidth A'-fmec'^. On inverse- Compton scattering the effective occupation number of 
scattered photons will be reduced by the ratio of bandwidths, (A'-f /'-/)/ {Au/u). If we take 
(A7/7) ~ 1, and (Ai//i/) ~ 3 x lO^'' corresponding to a 1 km s~^ velocity spread, the 
reduction in occupation number is of order 3 x 10^ bringing the actual occupation number 
down to ~ 10^. 

The occupation number of inverse- Compton scattered photons also could in principle 
be diluted by the low, effective cross section for back-scatter, i.e. by the Klein-Nishina cross 
section for back-scattering. However, despite the 7^ value of 10^(5/25) (7/(4 x 10^)) for 
electrons, the incident photons only have energy ^hu5 ~ 10^ eV in the electron's rest frame. 
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far lower than the 0.511 MeV electron rest mass. The Klein-Nishina cross section, therefore, 
reduces to the Thomson cross section gt ~ 6.6 x 10~^^cm^. 

We can assume that the masers are isotropic, or else, if they are not, that there are 
a larger number than are actually observed. Either way, the scattered light they produce 
would be the same. If we further assume a jet with relativistic electron column density 
through which the maser photons pass of order Ue^rei^ ~ 10^''cm~^, wc can estimate the 
phase cell occupation number of the scattered radiation. It is the product of the maser 
beam phase-cell occupation number, the ratio of bandwidths, the electron column density, 
and the Thomson cross section, giving 



^ \3x lO^V V 3 X 10-6 J\l0^^cm-y 

Interestingly, those phase cells with high back-scattered occupation number will 
increase their occupancy at a rate {Nh + 1) times faster than unoccupied cells, since induced 
scattering then begins to play a role — there is gain. We may, therefore, expect such a 
configuration to give rise to reasonably high occupation numbers for TeV photons and 
energy densities compatible with observed values. NGC 1052, exhibits nearly 40 maser hot 
spots, with a total 22 GHz luminosity of ~ 2OOL0 ~ 8 x 10^^ erg s"^. Let us assume that 
the maser power available for interacting with the relativistic jet would be equivalent to 
only 25% of this. If only fraction n,, ,.^,il(jT ~ 6.6 x 10"^ of this radiation is scattered, but 
each photon's energy increases by ~ 1.1 X IQi^ the 1 TeV luminosity is ~ 1.5 X 10 ergs . 
This needs to be compared to the TeV flux from Mkn 501 in its high state, which is 
of order 3 x lO"^'' ergcm"^ s~^, corresponding for a distance of 200 Mpc to an apparent 
omnidirectional luminosity of 1.5 x lO^^ergs"^ (Plan et al. 1998). Since our model assumes 
only a single jet spherically expanding within a relatively narrow solid cone whose axis is 
directed at us, these two figures are roughly consonant. 



7. Synchrotron Emission from Relativistic Electrons 

A highly relativistic electron with energy E emits synchrotron power in its rest frame 
P{E) ~ 2.6 X 10-^(5/0. lgauss)2(7/4 x Wfeigs'^. 

The peak frequency the photons attain in this frame will be of the order of 

v„~#^~4.5xlO-r^V^yHz (2) 



2'KmeC ' VO.lgauss/ V4 X 10^ 
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where m.e is the electron rest mass, B ~ 0.1 gauss (e.g. Bednarek & Protheroe 1999) is the 
local magnetic field strength, and e the electron charge. In the terrestrial observer's frame 
the frequency becomes ~ 10"^^ Hz, which roughly corresponds to the peak synchrotron 
radiation frequency of Mkn 501 in the high state. 

OSSE observations during flaring (Catanese et al. 1997) show that the energy flux per 
log energy interval continues up to ~ 500 keV at roughly the same level as that observed 
by Beppo-SAX (Plan et al. 1998), indicating that Mkn 501 emits a synchrotron power at 
0.5 McV comparable to the TeV power during flaring. The emitted synchrotron power in 
the rclativistic jet's comoving frame would be 2.4 x 10^^(25/5)^ ergs implying emission 
from ~ 10^^(25/5)2(0.1gauss/5)2(4 x lO^/jf rclativistic electrons. 

In recent models of AGN jet dynamics (e.g. Falcke & Biermann 1999) a rclativistic jet 
can readily interact with A^ci dense ambient molecular clumps located at ~ 10^^ cm from 
the central engine, to produce rclativistic shocks that could trigger maser emission in these 
clumps. Local acceleration at the shock fronts or production from hadronic interaction 
and decays could then also provide rclativistic particle energies jrUeC^ ~ 8.2(7/4 x 10^) erg 
in the jet's comoving system. The time scale for energy loss for these particles through 
synchrotron radiation is of order tgynch ~ 1-25 x 10^(4 x 10^/7)(0.1 gauss/S)^ seconds. Since 
the rclativistic shocks propagate into the jet at a significant fraction of the speed of light, 
the radiating post-shock volumes have dimensions of order ctgynch ~ 10^^ to 10^^ cm on a 
side. At particle densities of order Ue^rei ~ lO^/A^^d cm~^, a post-shock column density of 
10^^ cm~^, through the A^ci shocks, therefore, appears possible. 



8. Discussion 

It is possible that highly energetic gamma radiation from distant cosmological sources 
will be found to appear in conflict with pair-production constraints imposed by the diffuse 
extragalactic infrared background radiation. This apparent violation could then be due 
to coherent TcV gammas of lower energy, whose Chcrcnkov radiation superflcially mimics 
individual photons of much higher energy. We have suggested how the Cherenkov radiation 
signatures of coherent and incoherent radiation can be distinguished, and have sketched a 
plausible way in which coherent TeV photons could be astrophysically generated. Whether 
this particular mechanism is found in nature, remains to be determined, but other possible 
sources of coherent TeV gamma radiation are also entirely possible. If coherent TeV 
photons can be produced in nature then wc have shown that there exists a mechanism by 
which multi-TeV Cherenkov signals may be observed from high redshift sources. 
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Fig. 1. — Number of electrons, A'^e, in the air shower versus depth in the atmosphere measured 
in radiation lengths. The top three curves show Ne and ± la for single photons of energy 
25 TeV, the middle three curves show and Ne ± la for an assembly of five coherent 
photons each of energy 5 TeV (note these curves have been displaced down by a factor of 5 
to avoid confusion), and the bottom three curves correspond to an assembly of twenty- five 
coherent photons each of energy 1 TeV (displaced down by a factor of 25). 
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Fig. 2. — Propagation of coherent photons across 200 Mpc through the microwave and 
diffuse infrared extragalactic radiation fields. No is the number of photons in the coherent 
assembly emitted in our direction. Nji is the surviving number, ii^obs is the energy of the 
observed Cherenkov event, equivalent to NjiS, where e is the energy of individual photons 
in the assembly. 



